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Abstract 

We construct a left-right symmetric (LRS) model in five dimensions which 
accounts naturally for the lepton flavor parameters. The fifth dimension is 
described by an orbifold, S1/Z2 x Z2, with a typical size of order TeV^^. 
The fundamental scale is of order 25 TeV which implies that the gauge hi- 
erarchy problem is ameliorated. In addition the LRS breaking scale is of 
order few TeV which implies that interactions beyond those of the standard 
model are accessible to near future experiments. Leptons of different repre- 
sentations are localized around different orbifold fixed points. This explains, 
through the Arkani-Hamed-Schmaltz mechanism, the smallness of the tau 
mass compared to the electroweak breaking scale. An additional U(l) hor- 
izontal symmetry, broken by small parameters, yields the hierarchy in the 
charged lepton masses, strong suppression of the light neutrino masses and 
accounts for the mixing parameters. The model yields several unique pre- 
dictions. In particular, the branching ratio for the lepton flavor violating 
process ^~ e^e~e~ is comparable with its present experimental sensitiv- 
ity. 



I. INTRODUCTION 

The recent results from the SNO [1] and other [2] experiments provide strong evidences 
for the incompleteness of the standard model (SM). Among the various new physics (NP) 
scenarios that predict neutrino masses, the left-right symmetric (LRS) framework [3] is an 
attractive and popular one. 

In many LRS models, e.g. models embedded in GUT [4], the LRS breaking scale, 
vr ^ 10^'' GeV, is much higher than the electroweak (EW) breaking scale, ki ~ 10^ GeV, 
and the low energy effective Lagrangian is similar in many aspects to that of the SM. In 
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such a case present and near future experiments will not be able to directly probe the 
NP. Furthermore, the introduction of such a high scale, in addition to the Planck scale, 
Mp ~ 10^^ GeV, related to gravity in four dimensions, raises the gauge hierarchy and fine 
tuning problems shared by many models. 

A new exciting possibility, however, that the fundamental scale of gravity can in fact 
be much smaller than Mp was raised in [5] . It is very interesting, therefore, to investigate 
whether a natural LRS model (LRSM), in which both Mp and vr are low, say below 100 
TeV, and in which the neutrinos are very light, can be constructed. 

This work is focused on the lepton flavor parameters, we comment on the inclusion 
of the quark sector in the conclusion. We present a LRSM which naturally accounts for 
the flavor parameters of the lepton sector and in which the fundamental scale and the 
LRS breaking scale are of the order of or below 25 TeV. It is a model of, at least, ^ one 
extra compact dimension which copes with the above problems by using both the Proggatt- 
Nielsen [6] and the Arkani-Hamed-Schmaltz (AS) [7] mechanisms. The main role of the 
AS mechanism is to localize "left" and "right" lepton fields on different fixed points in the 
extra dimension, which then explains the smallness of rur compared with the electroweak 
breaking scale. An additional U(l) horizontal symmetry, described in detail in [8], yields 
a modified see-saw mechanism and accounts for the other fiavor parameters. 

Before proceeding with the details of our work we note that recent two papers, by 
Mimura and Nandi [9] and by Mohapatra and Pcrez-Lorenzana [10], dealt with the con- 
struction of a LRSM in five dimensions [5D]. Though important parts of the analysis 
in [9,10] are general and also used below the three models use, in fact, different construc- 
tions in the extra dimension. Furthermore, to the best of our knowledge, the present model 
is the only LRSM in 5D that aims to account naturally for all the lepton flavor parameters. 

To make our discussion concrete we briefly list the lepton flavor parameters as deduced 
from various experiments. The charged lepton masses are [11], 

me ~ 5.1 ■ 10"^ MeV, ~ 1.1 • 10^ MeV , ~ 1.8 • 10^ MeV . (1) 

As for the neutrino parameters, we consider below the large mixing angle solution of the 
solar neutrino problem which is favored by data [12,13]. Consequently, the neutrino mass 
differences, at the 3(7 CL, are: 

AmL = (0.2 - 3) • 10-^ eV^ , Ami,^ = (1 - 6) • 10"^ eV^ , (2) 

where Amgoi [Am^tm] the mass square difference deduced from the data of the solar 
[atmospheric] neutrino experiments. The neutrino mixing parameters are [12,13]: 



To have Mp ^100 TeV requires more extra dimensions. 
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tan^ 012 = 0.2 - 0.7 , tan^ 623 = 0.4 - 3 , ^13 < 0.15 



(3) 



In addition there are both direct [11,14,15] and indirect [16] bounds on the absolute scale 
of the neutrino masses: 

mi < 1 eV . (4) 

In section II we present our 5D LRSM model. In section III we describe the 4D effective 
theory of our model and calculate its lepton flavor parameters. In section IV we suggest 
several criteria by which we can test the model predictions. Comments and conclusions 
are given in section V. 

II. THE MODEL 

The space time of our model is described by the usual 4D space and an additional space 
dimension compactified on the orbifold S1/Z2 x The characteristic energy scales are: 



L-i~ TeV; L-'<vr,v<M,; a = v^—>5, (5) 

where L is the size of the fifth dimension fundamental domain, vr is the scale at which 
the LRS is broken spontaneously, v is related to the typical width of the fermion wave 
functions and is the fundamental scale. 

The discrete group, Z2 x Z2 yields the following identifications for the fifth dimension 
coordinate, y: 

Z2- y^-y; Z'^: y' ^ -y' , (6) 

where y' = y + 7iR/2. 

The symmetry of the model is given by, ^ 

SU(2), X SU(2)2 X U(1)b_l x Glr x U(1)h , (7) 

where Glr corresponds to parity symmetry in the usual 4D LRSM and U(1)h corresponds 
the global horizontal symmetry [8]. The gauge group is broken both explicitly, by the 
transformation laws of the fields under the orbifold discrete group, and spontaneously, by 
the VEVs of the scalars [9]. 



^To avoid confusion when dealing with fermions in 5D wc denote the two SU(2) gauge groups as 
SU(2)i^2- We switch back to the ordinary notations, S\J{2)l^r, when we discuss the 4D effective 
theory. 
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We now move to describe the field content of our model. We first describe the scalar 
sector, then we move to the lepton sector and afterwards to the gauge boson sector. For 
each sector we describe the transformation of the fields under the gauge and the orbifold 
groups while the horizontal charges are given in the next section. 

The scalar field content of the model is similar to the minimal LRSM (see e.g. [3,8,17]). 
01 is a bi-fundamental of the two SU(2) groups and Ai^2 are triplets of the SU(2)i_2 gauge 
groups: 

In addition, we assume the existence of a real scalar field </?, a pseudo-singlet of the discrete 
LRS group, 

(p ^ -(p^ (9) 

where its self interactions and coupling to the fermions are discussed in the appendix. 
The transformation laws of the scalars under the orbifold discrete group are given by, 

A2(x^ -y) = A2(a:^ y) , A2(a:^ -y') = -P'^^i^x^, y')P' ; 
Ai(x^-|/) = Ai(a:^i/), Ai(x^-y') = Ai(x^y') 
(^(x^ -y) - -(^(x^ y) , ^{x\ -y') = -(^(x^ y') , (10) 

with P' — diag(l,— 1). Note that, as already discussed in [9], only one of the neutral 
components of 0i, h\, has a zero mode and can develop a VEV. This fact is related to the 
natural suppression of the Dirac neutrino masses, which is necessary for the phenomeno- 
logical viability of our model [8] . 

The field content of the lepton sector is more involved. It is similar to the one of [10] 
but not identical since our mechanism of generating neutrino masses is very different from 
the one of [10]. In most of the LRS models there is a pair of lepton doublets (connected 
by Glr) for each generation. In our model, we actually introduce two such pairs for each 
generations (see also [10]). The reason for this is that our model assumes canonical see-saw 
mechanism [18]. This implies zero modes for both the 4D left and right handed neutrinos. 
The transformation law of the bidoublet 0i (10) requires that the doublets of the SU(2)2 
group will transform non-trivially under the orbifold Z2 discrete group. Thus, only one 
of the two components of the doublet can have a zero mode. This enforces to double the 
number of fermion fields as previously done in [10]. Consequently, for each generation i, 
we have four lepton doublets. Two doublets of SU(2)i and two of SU(2)2, denoted as L\, 
L'l and ^2) -^2 respectively. 
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The representation of the leptons under the SU(2)ixSU(2)2xU(1)b-l gauge group is 
therefore given by, 

Li(2,l)_i ^ L^2(l,2)-i, L'I{2,1U ^ m,^U, (11) 

where i = 1..3 stands for lepton flavors. 

The lepton fields have the following transformation laws under the discrete group: 

L{(x^, -y) = -Li(x^ y) , L{(x^, -y') = -U,{x\ y') ; 
L'iix^, -y) = -L^(x^ y) , L'/(x^ -y') = L'i{x>^, y') ; 

-y) = , 4(^^ -?/') = y') ; 

-y) = |/) , -|/') = -P'L'^{x^, y') , (12) 

where note that unlike [10] the transformation laws of L2 allow for the "right handed" 
neutrinos to have zero modes. As we shall see below, however, they acquire large masses 
due to their Yukawa couplings to A2. 

In order to have a chiral 4D low energy effective theory any lepton, ip, has, on top of 
(12), the following transformation law under the orbifold discrete group [19,20]: 

V'(x^ X,) ^ ij(x'', -X,) = ^'ij(x^, X,) , (13) 

where x^ stands both for y and y'. 

We now move to the gauge boson sector. As was already discussed in [9,10], with the 
above transformation laws for the matter fields, the gauge bosons of the SU(2)2xU(1)b_l 
gauge groups have the following transformations: 

W^{x>^, -y) = W^{x^, y) , W^{x^^, -y') = P'W^{x^, y')P' ; 
W,{x>^,-y) = -W,{x^^,y), W^{x^^,-y') = -P'W^{x>' ,y)P' ■ 
B,{x^, -y), B,{x^, -y') = B,{x^, y), B,{x^, y') ; 

B^{x^,-y),B^{x^,-y')^-B^[x^,y),-B^{x\y'), (14) 

where B^,^^ corresponds to the U(1)s_l gauge group. The gauge bosons of the SU(2)i have 
the following transformations: 

W^{x^, -X,) = W^ix^', X5) , W^ix^, -X,) = -W,{x^, X,) , (15) 

where X5 stands both for y and y'. 

The additional U(1)h symmetry (discussed in detail in [8]), is assumed to be broken by 
small parameters in two stages. First it is broken to a discrete Z2 symmetry (the discrete 
symmetry does not allow for Majorana masses) by a small parameter, e. Then the discrete 
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subgroup is further broken by a small parameter, 6. Thus, as discussed below, various 
terms in the 5D effective Lagrangian are suppressed by powers of e and 6: 

jC^D ^^^^ 1 J^|Q(4)-Q(Li)+Q(0OI /Zi0i4 + £l«(^'2')-Q(M)-Q{</>i)I^Zi(/.2Li' 

^ Se\Q(^^)+Q(LD+Q(^i)^h\h{Lfir2A,Li + Lfir2A2Li) + ^ 

+ V^'^(0,A), (16) 

where £o contains the kinetic and the gauge interaction terms [9,10]. Higher dimensional 
operators are subdominant for the low energy effective theory. This is due to suppression 
factors coming from inverse powers of ~ 25 TeV and from powers of e and 5 [21]. ^ 
Furthermore, we assume that quantum modifications to our model both from perturbative 
and non-perturbative sources are small in the IR hmit of the 4D effective theory (for 
discussions on this subjects see e.g [22,23] and references therein). 
The effective 5D scalar potential, V^^{(f),A) is given by 

V'^i^P, A) = ^ {-/.^,Tr(0l0,) - i^^Tr{AAl) + ^ [A,,Tr(0t0^.)Tr(0l0,) 

ijk * 

+ p,,Tr(A,A,)Tr(A]A]) + a,,Tr(0j0,)Tr(A,Al) + A,TV(A1 Ai0,0| 

+ A^A20l0,) + 72iTr(Al02A20l)] } , (17) 

where some of the coefficients in V^^ contain implicit suppression factor of various powers 
of e [8] . In the generic case, ^1.2 and Ai^2 develop VEVs, 







(18) 



The different VEVs have the following relation among them [3,17], 

\vrVl\ = ■ (19) 

In the presence of the U(1)h 7 is given by [8], 

7 g|Q(AK)-Q(Ai.)-2Q(</>i)| _ ^20) 



^As an example consider the operator L\L\(l)2/^2(ki which is similar to the operator that gives 
the dominant contribution for neutrino masses in [10]. In our model the induced neutrino mass 
from such an operator is roughly, e'^'^dklvji/M^, which is completely negligible. 
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III. THE SPECTRUM OF THE 4D THEORY 



Most of the details of the model were given above. To calculate its low energy spectrum, 
however, two main ingredients are missing. One is the horizontal charges of the fields, which 
are given below. The other is related to the separation and localization of the different 
lepton fields. We assume that all the doublets of the SU(2)i gauge group have same-sign 
Yukawa couplings to ip (A2). Since we choose ip to be odd under the LRS discrete group 
(9), all the SU(2)2 doublets have the opposite sign of Yukawa couplings to if. ^ 

In the model presented below, we assume that these Yukawa couplings are flavor in- 
dependent. Consequently, the magnitude of the couplings between the leptons and ip is 
universal and naturally of order unity. This assumption can be motivated in cases where 
the couplings to the scalar are originated from a flavor blind sector of a more fundamental 
theory. We shall comment on the implications of relaxing this assumption in the following 
section. 

As explained, the fields L\ and L\ are localized around different fixed points [19,24]. 
In the appendix we show that bidoublet Yukawa couplings arc naturally suppressed due 
to small overlap between the wave functions of the zero modes of L\ and As shown in 
(A21), given the model fundamental parameters (5), the corresponding suppression factor, 
K ~ e~", is naturally of the order of rrir/ki. 

At this stage, when our main focus is on the zero mode of the fields, it is convenient to 
switch to the ordinary 4D notations, SU{2)i^2 ^ SU(2)l^r; Ai,2 ^ ^l,r-, /i°2 ^ 
where /° is any fermion zero mode that belongs to a doublet of the SU(2)j gauge group. 

The charges of the fields under the U(1)h horizontal symmetry are given by: 

Q(Li) = o, g(Li) = i, g(Li) = 3, 

Q{Al) = -7/2 , Q{Ar) = -21/2 , g(0i) = 7 , (21) 

with Q{Ll) — Q{L'f). Note that the charges of Al,r are half integers and therefore they 
carry an odd parity under the residual, Z2, horizontal symmetry. All the other fields have 
integer charges and therefore carry an even parity. Thus, in the limit where the discrete 
symmetry is exact {S 0) the Majorana mass matrices vanish. The discrete symmetry is 
assumed to be further broken by the small parameter 6. 

For concreteness, in our calculation below, we use the following numerical values, 

£, 5 ~ 0.3 ; ~ 4 TeV . (22) 



■^In that way we go one step fiirther in eliminating the arbitrariness in the sign assignment of 
the corresponding Yukawas. Such a problem is often encountered in the AS framework. 
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Other combination of parameters of a similar magnitude yield a viable phenomenology as 
well. Furthermore, for simplicity we assume real values for the various VEVs and couplings. 
Consequently, we assume no CP violation in the lepton sector. 

We now arrive at a point where we can compute the masses of the various lepton zero 
modes predicted by the model. 



A. Charged leptons 



The charged lepton mass matrix, M'^, is read from the second term in the square 
brackets of eq. (16). It is given, up to order one coefficients, by: 



Kki 



(23) 



1 / 



Using the numerical value for Kki ~ 2 GeV, the eigenvalues of M^' reproduce the required 
scale for the charged lepton masses (1), up to order one coefficients: 



rUe ~ e^Kki ~ 1 MeV , ~ e'^Kki ~ 150 MeV , ~ X/ci ~ 2 GeV . 



(24) 



In addition M is hierarchical and diagonalized by. 



ei,r^e\ ei,^e\ ei,r^e. 



(25) 



B. Neutrinos 

The light neutrino mass matrix. Mi , is given by: 

~ {M^fiM^^^y^M^ + , (26) 

with being the Dirac neutrino mass and M^"^ [M^"^] being the Maj orana mass 
matrix for the right [left] handed neutrinos. The RHS of eq. (26) contains two terms. The 
first, M^see = (M^)^(M^"^)-iM^, is related to the seesaw mechanism [18]. The second, 
is induced by the VEV of A^. In the following we shall calculate the elements of 
each and we shall show that the dominant contributions to M/^ come from M^""^ , while 
M^gg yields small (but non-negligible) correction to Ml'. 

In our model M^"'^ is read from the = Ai Yukawa interactions of eq. (16). 
Since A^ r carry half integer charges their corresponding couplings are suppressed by 
5 X niax(£l«(^^.«)+2Q(^^.«)+V2|^£|Q(A^,«)+2Q(L^.«)-i/2|j_ Consequently, Mf"^' is given by. 
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e' 1 1 
M^''^ = vl5\ 1 e 



(27) 



1 £2 ^3 



where from eq. (20) we have, 

vl ~ 



M.=-IQ(Afl)-Q(Az,)-2Q(</>i)| _ M., 



.21 



(2^ 



vr Vr 

As explained above M^^g depends on the structure of the Dirac neutrino mass, M^f , 
and on the structure of the Majorana mass matrix for the right handed neutrinos M^"''' . 
The matrix M^""^ is simply read from M^""' with the replacement vl vr, 

£2 1 1\ 



vr5 



1 

1 £ 



e e 

2 ^3 



2 



(29) 



The matrix M^l^ has an approximate L^ — L^ — structure [25,26]. Thus to diagonahze 
^R,L requires 9x2 ~ 7r/4, ^23 ~ 1 and a very small ^13 (for a recent review see e.g. [13] 
and references therein). 

Since the determinant of M^"'^ is given by. 



Det {m^^') ^ {vR6 fe , 



the inverse of M^"'-' can be approximated by. 



(Mr^)-i 



VrS 



O23{023)Ou{ei2)diag{l, 1, 1/£)0^2(^12)0^3(^23) , 



(30) 



(31) 



where Oij{9ij) is a rotation matrix on the ij plane with an angle 9ij which were defined 
above. 

The neutrino Dirac mass matrix, M^, is read from the first term in the square brackets 
of eq. (16) . It is given by. 



(32) 







£^0 




~ Kki 












£l3 


£lV 



We can now use eqs. (31,32) to calculate M^g^, 

/I ^2 ^3' 

,5 fi-l 



1 £^ £=^\ 



S Vr 



£2 £4 £5 
,£=^ £5 £V 



1 


1 \ 


/I £2 £3 \ 


l/£ 




£2 £4 £5 









~ £ 



23 ""1 



1 £2 ^3- 



, £3 £5 £6 



(33) 
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where in the second hne we used the approximation, ~ 

Combining eqs. (27,33) the hght neutrinos mass matrix is given by. 



,2/^1 1 



1 e e'\ , (34) 
VI e 



where, as anticipated, the dominant elements come from M^"'-' . Consequently, the Matrix 
Mj" has an approximate L^ — L^ — Lj- structure which therefore yields inverted hierarchical 
masses: 

^2 ^2 
Imil ~ Imal 0.02 eV , \m^\ fa e^^— - 0.007 eV ; 

Vr Vr 

AmL ~ 2 • 10-^ eV^ , Am^^ ~ 1 • 10"=^ eV^ , (35) 

in agreement, up to an order one coefficients, with the recent data given in eq. (2). The 
neutrino mixing angles read from Mi and M'^^ (23) are: 



yi+tan2 2^12-1 
tan' ^12 = ' . =■ 1 - £ ~ 0.7, 

^1 + tan2 9^2 + 1 
tan^ ^23 = £>(!), ^13 ~ 0{e^) ~ 0.1 , (36) 

again in agreement with the recent neutrino data given in eq. (3). 



IV. TESTING THE MODEL BY EXPERIMENTS 

It is very interesting to understand whether the model has unique experimental signa- 
tures which can be tested in near future experiments. Our model belongs to a class of LRS 
models in 5D, that was recently constructed in [9,10] and some of its phenomenological 
implications were already discussed there. We shall briefly summarize some of the most 
important properties: 

(1) Existence of Kaluza-Klein (KK) excitations of the gauge bosons; (2) No left-right mix- 
ing in the charged sector; (3) The hghtest 1^2 is a KK excitation while zf^ is not; (4) A 
lower bound of the order of a TeV on and on the masses of the Z2 and 1^2 lower KK 
excitations; (5) Existence of a heavy stable lepton whose mass is of the order of L~^. 
In addition we list below predictions which are specifically related to our model: 

(i) Inverted mass hierarchy for the neutrinos. 

(ii) Order one mixing for ^23 but not parametrically close to maximal. 
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(iii) Small branching ratio (BR) of lepton flavor processes such as ^ 67, r ^ fi'-f and 
r ej (for details on these processes see e.g. [27] and references therein) due to 
the smallness of the LR mixing in the model [9,10]. The same feature is also shared 
by the additional neutral and charged scalars introduced in our model. This is due 
to their corresponding transformation under the orbifold discrete group (10) and the 
smallness of VL/{k\,VR) (20). 

(iv) The process /i~ — > e^e'e" does not require LR mixing and therefore might be signifi- 
cantly enhanced (similar r decay modes have a smaller experimental sensitivity [11]). 
It can be efficiently mediated by the Yukawa couplings of to a pair of charged 
lepton zero modes. The rate is given by (see e.g. [27-29] and references therein) 

where Ge// is proportional to the product of (^^'"ee) and {S2~eii) Yukawa couplings. 
In our model it is given by 

Comparing the above result with the current experimental bound [11], 

G'e//< 10-^1 GeV-2, (39) 

and we learn that the two are comparable. This implies that our model will be 
subject to an experimental test in the very near future. 

(v) As Ml (34) is induced by an approximate Le — L^j, — Lr symmetry, the model predicts 
a rather small amplitude for neutrinoless double decay (see e.g. [15,30] and references 
therein) : 

|(m)|-£2^ 0.007 eV. (40) 

This is comparable with the lower limit for the case of inverted hierarchy [15] and 
very likely, below the sensitivity of near future experiments. 

The above results are based on the assumption that the widths of the lepton wave 
functions are universal. If the Yukawa couplings between the leptons and (p are not universal 
but of order unity the model looses some of its predictive power. This is mainly due to 
the fact the structure of the Dirac mass matrices (23,32) may be significantly modified due 
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to different overlaps of tlie wave functions in tlie extra dimension (even if tlie eigenvalues 
of these matrices are unchanged). In this case the resultant mixing angles and neutrino 
masses are functions of the widths and cannot be explicitly calculated. Note, however, that 
if the hierarchical structure (23,32) is not badly broken (with the corresponding eigenvalues 
remain roughly the same), most of the above features of the model are still vahd. This is 
due to the fact that the structure of the Majorana mass matrices (27,29) is independent of 
the universality assumption and the neutrino masses and mixing are mainly derived from 
Mf (27). 

V. DISCUSSION AND CONCLUSION 

We constructed a left-right symmetric (LRS) model of leptons in five dimensions. The 
structure of the lepton flavor parameters is explained naturally and the gauge hierarchy 
problem is ameliorated. 

The typical scales are: (1) The fundamental-domain size in the extra dimension, L ~ 
TeV"^; (2) The LRS breaking scale, vr A TeV; (3) The fundamental scale, ~ 
25 TeV. 

The model predicts that the light neutrino mass matrix has an approximate Lg — L^ — L,- 
structure which therefore yields inverted hierarchical masses, with tan^ O12 ~ 0.7 and ^23 — 
0{1) but not parametrically close to maximal. The model predicts a relatively large 
branching ratio for the lepton flavor violation process — > e^e~e~ , comparable with its 
present experimental sensitivity. 

The fact that the model copes both with the lepton flavor puzzle and with the gauge 
hierarchy problem is manifested in its rich structure. The model has elements related to 
the following mechanisms: 

(i) The Arkani-Hamed-Schmaltz mechanism- localization of "left" and "right" leptons on 
different fixed points in the extra dimension, which then explains the smallness of mi- 
compared with the electroweak breaking scale. 

(ii) The Froggatt-Nielsen mechanism- a horizontal symmetry which yields a modified see- 
saw mechanism and accounts for the other flavor parameters. 

The above results are derived for the case of universal Yukawa couplings to ip and when 
L* and L'* have the same horizontal charges. If this is not the case the model looses some 
of its predictive power. It is rather clear, however, that it is possible, for a given set of 
order one Yukawa couplings, to find a charge assignment that preserve the model appealing 
features. Thus, although the more generic setup of the model is less predictive, it is still 
provide a natural framework to understand the origin of the fiavor parameters. 

In this work we focused on the large mixing angle solution of the solar neutrino problem 
as favored by the present experimental data. In an earlier work [8], however, we showed 
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that a LRSM with a similar and the same horizontal symmetry, can in principle, produce 
a model with the small mixing angle solution of the solar neutrino problem. 

Finally we note that the addition of the quark fields to the fermion sector of our model 
is in principle possible based on the above ideas [24]. The large top mass, however, might 
require an additional structure [31]. 
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APPENDIX A: LOCALIZATION OF THE FERMION WAVE FUNCTIONS 



In this part we examine the requirements for sphtting and locahzing the fermions around 
the different fixed points. Some of the results derived below are known (see e.g. [7,19,32]), 
but for the self consistency of our work we briefly present them below. 

To localize the fermions we add a real scalar field, </?, which transforms non-trivially 
under the orbifold discrete group [7,19], 

ip{x,X5) ^ ip{x,-X5) = -ip{x,X5) , (Al) 

which for simplicity we take as a single Z2. The Lagrangian is 

£ = V - 75^5 - V' + Id^^d^^ - \d,vd^V -X{^^- v'^)" , (A2) 

where for further simplicity we consider a single fermion model. 

The shape of the fermion wave functions depends on the shape of {(p). The localization 
is significant if / dx^{(p{x^)) is large enough. In the following we investigate what is required 
from the model parameters so that efficient localization will occur. 



1. Properties of the profile of {ip{x^)) 



The boundary conditions, (Al), require that the scalar field vanish on the orbifold fixed 
points, namely, at = and L. However, for v''^ > 0, the scalar field tends to develop a 
VEV. We can make this argument more quantitative by considering the scalar Lagrangian, 
neglecting the effects of its Yukawa interactions: 



^d^(pdi^(p - ^d5(pd5(p - A((^^ - V 



I2\2 



(A3) 



The minimum energy configuration is found when the following action is minimized, 

.2 .1 n 



S=^J du ^d,^'d,^' + a\^'' - If 
L'^ Jo L2 



(A4) 



with u = x^l L, V = w'a/L, A' = AM^, a? = X'vL^ and Lp' = ^ip. ^ Let us find the 
condition for which p> = will not be a local minimum of S. Expanding (p in terms of sin 
functions which satisfy the boundary condition of eq. (Al), 



^Note that above (5) we defined = vL^. The difference is the A' coefficient which is assumed 
to be of order one and therefore it is omitted below. 
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and integrating over u yields the resultant "mass" matrix for the action 5", 



nm 



' 2 2 



Prom the structure of we learn that if 



(x6nm[^^-2a'\ . (A6) 
^(=0 V ^ / 



~ vL— > — (A7) 

we find negative eigenvalues for which implies that </?' = is not a local minimum of 
5" [19]. Thus, we expect to find a lower action 5" for u dependent (p', 

{^{u)) = h{u) . (A8) 

An intuition on the possible profile of h{u) is gained by viewing the action S as related 
to a one particle problem. The particle moves in a potential V{ip') = — a^(</?'^ — 1)^, 
shown in figure 1, with "initial" and "final" conditions y:>'{u = 0, 1) = (for more details 
sec e.g. [32]). The solution (p'{u) = represents a static particle at the bottom of the 
potential. 

We prove below that the "time" for the particle to reach the turning point and returning 
is a monotonically increasing function with the "length" of the corresponding "trajectory". 
Consequently, there is only one additional solution which must correspond to the global 
minimum, provided that eq. (A7) is satisfied. 

Since S does not contain an explicit u dependence, we can define an effective conserved 
"energy" , E, 

E=^d,^'d,^'-a\^"-l)\ (A9) 

where — < £^ < 0. Using this we can estimate the maximal value of (p' (the value at the 
"turning point" ) and study its dependence on the time period, T: 

f 1 fp'im d(p' 

1/2 



4 V2a Jo ^E/a? - {(p'^ - If 

1 ffi'im dip' 



V2aJo ^[cp'^-cp'{l/2n<p'^-^'{l/2y] 



= 1, (AlO) 



with v?'(l/2)^ = 1 — \J —E/a?, is the physical turning point and (^'(1/2)^ = 1 + \j—Eja^ — 
2 — (^'(1/2)^ is an unphysical one. 
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Using the substitution X = (^'/(^'(l/2), / takes the form of a complete ehptic integral 
(see e.g. [33]), 

dX _ F(n/2,d) ^^^^^ 



sj2o?ip'{l/2y Jo yJ[X^-l][X^-l/d''] sj2o?(p'{l/2y ' 

with d"^ = (/9'(l/2)^/(^'(l/2)^ and F{7r/2,d) is a complete elliptic integral of the first kind. 

The function F can be expressed as a series in d"^, (provided that < li < 1 as in our 
case) 



F(n/2,d) r^- {1+ , , , 
^ ' ' ' 2 1 Vl! -21 




(2n-l)!!^2n 
n!2" 



2" 

d'' 



(A12) 



As a consistency check for our results we note that for a very small (p' the potential is 
quadratic in (p', V{(p') ~ 2a^<^'^. This implies that the time period should be given by 

f^n/a. (A13) 

Plugging the first term in the expansion of F(7r/2, d) into the expression for T in eq. (AlO) 
indeed reproduces the result of eq. (A13). 

We are now at a point were we can show that the "time" for the particle to reach the 
turning point and returning is a monotonically increasing function with the "length" of the 
corresponding "trajectory" or "energy" . Since all the coefficients of the series in eq. (A12) 
are positive and d is monotonic with E (in the physical range) it is clear that F{7i/2, d) is a 
monotonic increasing function with E. The time period is given by / oc F{7i/2, d)/ip'{l/2). 
Since ^'{1/2) decreases with E, I or the time period is indeed increasing with the energy 
which guarantees that there is only one additional minimum of S. 

For the above potential, V{ip), one cannot compute analytically the profile of {ip{x5)). 
However to our consideration we only need to have information on its asymptotical be- 
havior. This was analyzed in [19] in which it was shown that for large L the profile is 
approximated by, in our scaling: 

h{u) a; tanh (^\/2auj tanh \/2a(l-'u) , (A14) 

where the exact u dependence of (A14) is found by taking the large a limit of eq. (A9). With 
that profile we expect a significant reduction in the overlap between the wave functions 
localized at the different fixed points. 

We demonstrate that this case corresponds to the large a case using our semi-analytic 
derivation presented above. In addition we gain physical understanding for the resultant 
shape of h{u). Consider the result of eqs. (AlO, All) when d --^ 1 or (/)'(l/2)^ ~ 1 — In 
that case F{t:/2, d) can be approximated by the following expression [33]: 
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F(7r/2,ci) ~ln4-lnV2£. 



(A15) 



Substituting this into eq. (AlO) yields the following equation for ip'(l/2) : 

a/V2 = ln4-lnV2£, (A16) 

which implies 

(^'(1/2)2 = l-Se-'^^. (A17) 

Prom (A17) we learn that indeed this hmit is good for a relatively large a. This means 
that at ~ the kinetic energy is very large and there is a very rapid motion there. 
Furthermore, since near the turning point the potential is almost fiat the net force on 
the particle is very small. This implies that the particle spends most of its time near the 
turning point moving with a nearly constant velocity. Thus, as promised, h{u) can be well 
approximated by the two kink profile given in eq. (A14). We also learn from (A17) that 
this limit is good for rather modest values of a. For example a > 4 yields 9?'(l/2)^ > 0.97. 



2. The shape of the fermion wave functions 



We analyzed above the profile of ((/^(xs)), relevant for our considerations. In this part 
we ask what are the conditions for the fermions to be localized around the different fixed 
points with small overlap of the corresponding wave functions. Using the derivation of 
[7,19] we find that the shape of a zero mode fermions in the fifth dimension is given by: 

eo+(x5) = K^,^e-<^'^ (A18) 

with K^Q^ being a normalization constant and 

six^ — f . — / duhiu) — fa I duhiu) , (A19) 



M 

where, as shown in [19], if / is positive [negative] the corresponding fermion will be localized 
around X5 = [L\. As explained above the fact that </? is odd under the LRS discrete group 
(9) guarantees that the fields Li, L\ [L2, L'^ have positive [negative] Yukawa coupling /. 

In order to estimate the overlap between the wave functions we approximate h{u) by 
a step function [24], h{u) ~ 1, in the required range. Consequently, the wave functions of 
two fermions with opposite Yukawas are roughly given by 

^1,2 ~ iri,2e^l-^i'^l"" . (A20) 

Thus, the overlap between the wave functions is given, up to an order one coefficients, 
by [24]: 
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K= ^7/^1-02 ~ e-^" , (A21) 
Jo 

where C is an order one coefficient, that for simphcity, we assume to be one. In the range 
a = 5.. 9 we get that the overlap between the wave functions is indeed small, as required 
by our model, of the order of rrir/ki. 
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FIGURES 
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FIG. 1. The black line corresponds to the effective potential which determines h{x^). The 
horizontal gray line represents a typical "energy" of h{x^)^ the "trajectory" which minimizes the 



effective action S of eq. (A4). The gray point, at 
solution. 



[0, —1], corresponds to the "static" 
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